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Abstract 

In a recent study (Jain et al 2007 Phys. Rev. Lett. 99 190601), a symmetric exclusion process 
with time-dependent hopping rates was introduced. Using simulations and a perturbation theory, 
it was shown that if the hopping rates at two neighboring sites of a closed ring vary periodically 
in time and have a relative phase difference, there is a net DC current which decreases inversely 
with the system size. In this work, we simplify and generalize our earlier treatment. We study a 
model where hopping rates at all sites vary periodically in time, and show that for certain choices 
of relative phases, a DC current of order unity can be obtained. Our results are obtained using 
a perturbation theory in the amplitude of the time-dependent part of the hopping rate. We also 
present results obtained in a sudden approximation that assumes large modulation frequency. 
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I. INTRODUCTION 



The symmetric exclusion process (SEP) is one of the simplest and well studied models of a 
stochastic interacting particle system. In this model which can be defined on a (i-dimensional 
hypercubic lattice, particles move diffusively while satisfying the hard core constraint that 
two particles cannot be on the same site. A number of exact results have been obtained for 
this model, particularly in one dimension If the model is defined on a ring and 

conserves the total density, the system obeys the equilibrium condition of detailed balance 
in the steady state and thus does not support any net current. A lot of attention has also 
been given to non-equilibrium steady states of driven SEP in which the particles can enter or 
eave the bulk at the boundaries. For this model, the time-dependent correlation functions 
J] and dynamical exponents have been obtained using the equivalence of the transition 
matrix (ly-matrix) to the Heisenberg model Recently, large deviation functional and 
current fluctuations have also been calculated for the driven SEP jo, 7, 8]. Experimentally it 
has been shown that SEP can be used to model the diffusion of colloidal particles in narrow 
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shown that similar effect can occur in a SEP model in which the hopping rates at two 
neig hbouring sites are chosen to vary periodically in time and with a relative phase difference 
291]. Our results obtained using Monte-Carlo simulations and a second-order perturbative 
calculation in the amplitude /i of the time-dependent part of the hopping rate can be 
summarized as follows: (i) A DC current J is obtained, which decays with system size L as 
J ~ 1/L. Correspondingly the time averaged density profile varies linearly in the bulk of the 
system, (ii) The DC current J depends sinusoidally on the phase difference between rates 
at two sites, (iii) The dependence of J on driving frequency u shows a peak at a frequency 
uj* with J^l/u;asu;^oo and J ^ u; as u; ^ 0. The latter result means that a finite 
number of particles are circulated even in the adiabatic limit. 
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es and heat, in similar time-dependent stochastic models, has 



331 ] and seen in experiments 3J]. Systems exhibiting pump- 



ing effect have often been modeled as Brownian ratchets in which non-interacting particles 
move in an external periodic potential. As discussed in 35|, these pump models are similar 
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to Brownian ratchets 



36l ] where non-interacting particles placed in spatially asymmetric 



potentials that vary periodically in time and acted upon by noise execute directed motion. 



Models o 
ered 37, 



' non-interacting particles moving in symmetric potentials have also been consid- 
and pumping demonstrated. However for the model studied by us, particle 



interactions seem necessary for the pumping effect . Our model differs from such models in 
that here we are dealing with a many body particle system with interactions. For such an 
extended system, as described in the following section, the n-point equation does not close 
and involves next order correlation functions also. Pumping effect has also been found in the 
steady state of a driven SEP with two species {A and B) of particles in which both species 



have the same diffusion constant 



391 ]. In this case, although the total current J a + Jb due 



to both species obeys the Pick's law, the current due to one of the species does not follow 
the density gradient. However the pumping mechanism is different from that in our model 
where it arises due to the time-dependent rates. 

In this paper, we consider a generalization of our earlier model by allowing the rates at all 
the sites to be time-dependent with a relative phase difference between neighbouring sites. 
The model is treated analytically using two approximations: a perturbation theory in the 
time-dependent part of the driving, and an expansion in the large frequency limit to leading 
orders. The treatment in this paper considerably simplifies the earlier one given in 29|. The 
most interesting new result is that in the model with time-dependent rates at all sites, a 
current of order unity can be obtained even in the thermodynamic limit for certain choices 
of relative phase differences. 

The paper is organized as follows. In section ((TTl), the model is defined. In section ( 1IIII) the 
details of the first perturbation theory (expansion in /i) are given, and two special choices of 
hopping rates are discussed. The results obtained from a sudden approximation (expansion 
in 1/uj)) are given in section (11 Vp . Pinally we end with a discussion in section (IVll . 



II. DEFINITION OF MODEL 



The model is defined on a ring with L sites. A site / = 1,2,3, ...L can be occupied by 
ni = or 1 particle and the system contains a total oi N = pL particles where p is the total 
density. A particle at site / hops to an empty site either on the left or right with equal rates 
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given by: 



ui = fo + flVl 

where vi = ai sm{ujt + (pi) = uie'"^* + i/fe"^'^* . (1) 

Here the site-dependent complex amphtudes are defined by ui = ai e*'^'/2z with ai real and 
/i is chosen such that all hopping rates are positive. We will discuss two particular choices 
for the hopping rates in detail. Our first choice corresponds to the case where the hopping 
rates are time-dependent at only two sites of the ring, and we get an average current which 
decays inversely with system size. In the second case, we choose time-dependent hopping 
rates at all sites and show that a finite current can be obtained even in the thermodynamic 
limit. 

A configuration of the system can be specified by the set {ni}, / = 1, 2, ...L. Let us define 
P{t) as the probability vector in the configuration space, with elements P{C,t) giving the 
probability of the system being in the configuration C = {rii} at time t. Then the stochastic 
dynamics of the many particle system is described by the master equation: 

= W{t) P(t) = Wo P(t) + Wi(t) P(t) (2) 

where W is the transition matrix, which we have split into a time- independent and a time- 
dependent part. One can also consider the time-evolution equations for m-point equal- 
time correlation functions Gi,i2,/3,....,i^(t) = {ni^...ni,J = Y.^ni}^h---^i,nP{{^i}^t)- Thus, for 
example, the density pi{t) = (rii) and the two-point correlation function Ci^mif) satisfy the 
following equations: 

dpi 
dt 

dCu 



+ 2niPi - ui_ipi_i - ui+ipi+i = ui{Ci_i^i + Ci^i+i) - ui+iCi^i+i - (3) 



+ 2{ Ui + Um )Ci m — Ul^i Ci^i m — Q+l,m ~ ^m-l C*/,m-l — ^m.+l Q,m,+1 

at 

ui ( Ci + Ci 

,l+l,m 

) +Um {Ci 

,m—l,m 

+ Ci Ci -l,l,m — 

l,m 

^■^m-l C«,m-l,m ~ '^m+l C*«,m,m+l5 ^OT \l — m\ 1 

+ {Ui + Ui+i - Ui^i - Ui+2 Ci^i+2 

Ul Cl-i^l^m + Ul+i Cl^l+i^l+2 — Ul-i Cl-i^l^l+i — Ul+2 C/,/+l,/+2 • (4) 



From Floquet's theorem 40|], it follows that the long time state of the system (assumed 



to be unique) will be periodic in time, with period T = 2-71/00. Here we will be mainly 



interested in the DC current J defined as 




(5) 



J Jo 

where the current Ji,i+i in a bond connecting sites / and Z + 1 is given by 



Ji,i+i — ui{pi — Ci^i+i) — ui^i{pi^i — Ci^i+i) 



(6) 



and the local density pi = (ni). From the periodicity of the state at long times and particle 
conservation, it follows that the DC current is uniform in space and therefore, using Eq. (I6l), 
we can write for the DC current: 



Thus to find the DC current, we need to compute the two-point correlation function C;^/+i(t). 
In this paper, we will develop two different perturbation schemes, valid for general vi, and 
then apply them to some special choices of the rates vi. 

Note that for /i = 0, the above model reduces to the homogeneous SEP with periodic 
boundary conditions whose properties are known exactly. In this case the steady state is an 
equilibrium state which obeys detailed balance and hence the average current is zero (note 
that this result holds even when the rates {ui} are site-dependent, but time-independent). 
In the steady state, all configurations are equally probable i.e. P{C) = 1/ (^) when /i = 0. 
Then one can show that the density and correlation functions for the homogeneous SEP are 
given by: 




(7) 





Pi 



,(0) 



P 



N 

T 




p 



(L-1) 




(9) 



III. PERTURBATION THEORY: EXPANSION IN / 



For fi 7^ 0, the knowledge of the exact steady state of homogeneous SEP enables us to set 
up a perturbation expansion in /i of various observables. We now describe this perturbation 
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k of the system. 



411 ] ■ We expand 



theory within which we calculate an expression for DC current J in the bu 
A similar perturbation technique was developed for a two-state system in 
various quantities of interest with fi as the perturbation parameter about the homogeneous 
steady state corresponding to fi = 0. Thus we write 

oo 

p,it) = {n,{t))=p + Y,flpt\t) (10) 

r=l 

oo 

cut) = {ni{t)nUt)) = C + E ' (11) 

r=l 

and similar expressions for higher correlations. Plugging in Eq. fill I) into Eq. ([H]), we find 
that the lowest order contribution to J is at 0{ff) and given by: 

j^'^ =7^1 - ^^+1 ) ^mIi • (12) 



1=1 



To develop our perturbation theory and find two-point correlation function C;*"^, we start 
with the time evolution equation for density pi{t) which is given by Eq. Q. Plugging in 
the expansions in Eqs. (fTU]) and (ITTl) . we get the following equation for the density pf'' at r^^ 
order: 



dp\ 



JO^lPl + mPi - Vl-^lPl^i - Vl+lPi^i 



dt 

= Mctif + ct;;^) - vi^ic^^^f - vi^ici:-i\ (13) 

where AiQi = gi^i — 2gi + gi_i defines the discrete Laplacian operator. Thus the density at 
^th Qj-f^QY f^Q^j^ obtained in terms of density and two point correlation function at (r — 1)*'* 
order. We check that at the zeroth order, we obtain the homogeneous SEP for which the 
density and all equal time correlations are given by Eq. (Q. At first order, the above equation 
then gives: 

"Pi f A .(1) 



at 



foAip\'' = roAiVi, (14) 



where tq = p — Cf'^. The solution for this equation is the sum of a homogeneous part which 
depends on initial conditions and a particular integral. At long times the homogeneous part 
vanishes while the particular integral has the following asymptotic form: 

pf\t) = Af^'"^' + Af^-'^^K (15) 
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Substituting Eq. (fT5l) in Eq. (fT^ we obtain the following equation for {A^^-*}: 

(zu + 2/o)A;^) - /oA[i\ - /o<\ = ro(z/,+i - 2z/; + z/,_i) . (16) 
This can be written in matrix form as: 

Z{uj) A = -ro B (17) 

where 

Zlm = —fo Sl,m+1 + { iuJ + 2/0 ) 6i^rn " /o <^/,m-l 

A = ^ = {^1, ^2, ^^l}^, (18) 

and periodic boundary conditions are implicitly taken. The above equation can be solved 
for A and we get: 

A = -To G{uj) B (19) 

where G{uj) = Z~^{uj). Both G{uj) and J3 are cyclic matrices and so can be diagonal- 
ized simultaneously. The eigenvalues of Z{uj) are iu + 4/o sin^(p7r/L), while that of B are 
4 sin^(j97r/L) with p = 1,2, ...,L, and eigenvector elements are q^'^'^p^/^ j i}/"^ _ Hence Ap"* can 
be written as: 

,(1) _ V- V- e L sm (pTT/L) 

A - L zo; + 4/osin2(p7r/L) ™' ^ 

m=l p=l \i / / 

which in the large L limit gives: 

A« = + ^ x: [ -L"-" + z^'-'i ] .m, (21) 

/o /o 2;+ - 

where, 2_ = y/2 - [{y/2f - lY'^ z+ = and y = 2 + {lu/U). 

To compute the 0{fl) contribution to J, we need to evaluate C/^, which we now proceed 
to obtain. Inserting the perturbation series in Eqs. ffTOj) and f|TT]) into Eq. (jl]) we get the 
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following equation for the correlation Ci^^ at r*'* order for |m — / | 7^ 1: 



dt 

- 1), ( r^''-'^'^ -L r'^'""^) ] 4-1, ( r'('"~^^ -^ r^'""^) 

Mr-l) Mr~l) r^ir-l) r'ir-l) 

while for m = / + 1 : 

M+i '-'i-i.i+i '-'«,«+2 



at 



Vl+2 \ ^1,1+2 " ^1,1+1,1+2 ) + I '-'1-1,1+1 " '-'1-1,1,1+1 I 

vi ( Cl^i+i - cl_^Jl_^_^ ) - vi+i ( Ci^i^j^ - Cl^i+i\+2 )• (22) 



At first order we get: 



+ /o (2C« , - C[%^, - C[%,) = Uv,^, + vi^, -V,- vi^,\ (23) 

where ko = C'/^/^ ~ ^^'^I'L '3 these are known from Eq. (Q. The computation of even the 
homogeneous solution of the above set of equations is in general a non-trivial task because 
of the form of the equations involving nearest neighbor indices and requires a Bethe ansatz 
or dynamic product ansatz 3, ^. However it turns out that the long time solution can still 
be found exactly and is given by: 

Ciiit) = ^fy^\t) + p^^\t)] = 4) + <^e-- , (24) 

where A^^^i = (^o/''"o)(A^^^ + A^). It is easily verified that this satisfies Eq. ( l23l) for all /,m. 
To determine whether the system indeed has a product measure requires a more detailed 
analysis of the higher order terms in the perturbation series and higher correlations. We 
have verified that at least to first order in perturbation theory, all correlation functions in 
fact have the same structure as the two-point correlation function in Eq. (124|) . 

We now plug the solution in Eq. flMl) into Eq. f|T2|) for the average current in the system 
and after some simplifications obtain: 

J'^' = ^ ( <S + <i-r - Af^'^^.i - A^^U, ) , (25) 
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FIG. 1: Log-log plot of J versus system size L at half filling for two cases discussed in the text. 
The points (♦) correspond to simulations for the 2 special sites problem with (p = 7r/2, /o = 0.3, 
/i = 0.2, u) = 0.27r. In this case the current goes as ^ . The points (•) corresponds to all 
sites having time-dependent hopping rates with q = -k 12 and /o = 0.5, f\ = 0.1, uj = 0.27r. In this 
case the current goes to a constant value at large L. The bold lines indicate the results from the 
perturbation theory. 

with a[^^ given by Eq. ( !2T]) . For any given choice of the rates ui, this general expression 
can be used to explicitly evaluate the net DC current in the system. We now consider two 
special choices of the rates {i^i}. 

(i) The choice ai = = 1, all other ai = 0, and 0i = 0,0^ = corresponds to the 



pumping problem with two special sites studied in 



29|. In the limit of large L, this gives: 



J(2) ^ ('A 
,/0 



kouj sin I 



Re[z_ 



(26) 



which agrees with the result presented in 29|] ( apart from a factor of two which was missed in 
that paper). Writing 2;+ = re*'', we find that for uj <^ uj* = 2/o, the magnitude r ^ l + y^uj/uj* 
and the angle rj ^ 



uj/uj*. In the opposite limit, r ^ lujjuj* and 77 ^ tx/2 — uj* /uj. Using 
2;+ = l/z_, we find that the current has the scaling form: 



J(2) 



flh sin i 
foL 



(27) 
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FIG. 2: Plot of DC current J versus density p = N/L for parameters /o = 0.5, /i = 0.4, lo = 0.27r 
and (pi = ttI/2 for system sizes L = 16,32 and 64. Both the results from simulations (symbols 
connected by dotted lines) and from the perturbation theory (lines) are plotted. 



where the scaling function G{x) = 2x for a; <^ 1 and 1/x for x ^ 1. We note that J is 
independent of /o for large x. This can be seen by writing the master equation as: 



dP 



— WoP(t) + — WiP(t) 



(28) 



d{ujt) UJ 

For UJ ^ fi), the first term on the right hand side can be neglected thus giving the probability 
distribution to be a function of fi/uj. 

(ii) The second case we consider here assumes a; = 1 at all sites and 0/ = g/, where 
q = 2tts/L with s = l,2...L/2, so that there is a constant phase difference q between 
successive sites. In this case, A\^^^s given by Eq. (120!) . evaluated at large L gives: 



where a 



2/o 
1 — cos q 

y/2 — cos q 



(29) 
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and from Eq. (l25|l we get for the average current: 

J*-^-* = — ^^j, ° sing Im[a] 
Jo 



2 /f /co sing (1 - cosg) ^^^^ 



2 ^ 4/2 (1 -cosg)2 



Thus we see that for most values of g we get a finite current, even in the hmit L — oo. For 
g ~ 1/L and g ~ tt — 1/L, the current goes to zero for large system size as J ~ L~^. From 
the current expression in Eq. (130|) . we can find out the value q = q*, at which the current is 
a maximum. By differentiating Eq. fl30|) with respect to g we get: 



cos(g*) = (1 + n^) - ^/ii + n^y-{i-n^), (3i) 

where Q = c<j/2/o. It turns out that for large u the maximum is at g* = 27r/3, while for 
small frequencies we get g* ~ ^/cJ. Also we find from Eq. (|30l) that in the adiabatic and fast 
drive limits, the currents are respectively given by: 

j(2)^|#^°*^^/2)^ ' ^//o«(l-cosg) ^^^^ 
2/2 /co sin g(l - cosg) ^ , ^//o » 1 . 

The perturbation theory results described above turn out to be quite accurate as can 
be seen from the comparisons with simulation results shown in Fig. ([1]) for both cases (i) 
and (ii). In this figure, we have plotted the current for different system sizes and verify the 
J ~ dependence for case (i) and J for case (ii) with g = 7r/2. Using the expression 
for /co in Eqs. (123 I5U]) . we find that J*^2) ^ p2^]^ _ which has a maximum at p* = 2/3 
and breaks particle-hole symmetry. This particle-hole asymmetry can be understood easily. 
From the definition of the model we see that, unlike the particles, the hopping rates of a 
hole are not symmetric: a hole at site / hops towards right with rate Ui^i and left with 
ui^i. In Fig. (I2l) we have plotted simulation results for the average current as a function of 
particle density, for different system sizes, and find good agreement with our perturbative 
result, even at a relatively large value of /i//o- 

In simulations we have looked at the density profiles and find that the site wise density 



profile pi in case (ii) is fiat. This is unlike in 29|], where we found high densities at the two 



special sites and then a linear density profile in the bulk. The fiat density profile, for case 
(ii), is understood because here there are no special pumping sites. It is interesting that we 
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can get current in the system even in the absence of Pick's law. We also note that even if 
the hop-out rates are made biased in one direction, like in the asymmetric exclusion process 
(ASEP), we can still get a current opposing this bias (for small biases). 



IV. SUDDEN APPROXIMATION: u;//o > 1 

In this section, we find the DC current within sudden approximation following the pro- 
cedure of 42|. Calling 9 = ut, the master equation Eq. ([2]) can be rewritten as 

^ = l[Wo + Wi(^)]P(e) (33) 

au to 

which can be expanded in powers of 1/uj by using P{0) = '^'^^q oo~^Pi"'\6) to give 

^?|^-Wi(»)Pf = W„Pi»> (35) 

au 

and so on. Prom the zeroth order equation, we see that Pf'^ is independent of 6. In fact, for 
u; — i> oo, we expect the system to behave as the unperturbed homogeneous SEP for which 
WoPi°'' = is satisfied and as discussed in Section [TTl all the elements of the vector P^"^ are 
known. Using this fact, the first order correction Pi^'' can be found by integrating Eq. 
over 6. Pollowing steps as those leading to Eq. f[T^ . we can obtain an expression for average 
current Js at order 0{l/u!) which is given by: 

E("« - (36) 

-^•J 1=1 

where we have expanded the nearest neighbor correlation function Ci^+i = Xl^^o ^~"'^u+i 
in powers of 1/u; and used the expression for Clf_^_i = C/^^^ given by Eq. ([9]). The first order 
correction to correlation function can be obtained by perturbatively expanding Eq. (j4j) and 
obeys the following simple equation: 

= fiko {vi+2 + vi^i- vi - vi+i) . (37) 
We now again discuss the two special choices of rates vi, discussed in the previous section. 
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(i) In this case, only two sites have time-dependent hopping rates. Solving the equations 
above for the correlation function, we get: 

C[^l = fMcos{9) - cos{9 + 0)) + ci,2 (38) 
C'i'li,^ = -/iA;o(cos(^) - cos(^ + 0)) + CL-i,L (39) 

c^l = fM^osie) + cos{e + 0)) + CL,i (40) 

where Cjj are constants of integration (which do not contribute to current). Using the above 
equations in the expression for Js^\ we finally obtain 

j(i) = 2/ffco^sin0 ^^^^ 

Thus, we find that to leading order in 1/ci; (and arbitrary /i), the DC current is the same 
as the one obtained by taking large lj limit in the current expression Eq. (1271) obtained from 
the /i expansion. 

(ii) In this case with a; = 1 at all sites, the equations for the first order correlation functions 
can be solved for arbitrary phases (pi, and we get: 

Cl}l, = koh [cos(^ + 00 + cos(e + 0i+i) - cos(^ + 0i_i) cos{e + 0^+2)] . (42) 

Using these in the current expression and after some simplifications, we get: 

= Ml ^ [2 sin(0,+i - 00 - sin(0,+i - 0,„i)] • (43) 



1=1 

Note that the above expression depends on the phase difference between nearest and next 
nearest neighbor sites. For 0; = g/, we recover the result stated in the second line of Eq. fl32l) . 



V. DISCUSSION 



In this article, we have considered a lattice model of diffusing particles with hard core 
interactions and shown that if the hopping rates at various sites are chosen to be symmetric 
but time-dependent, a DC current can be generated in the system. Thus a ratchet effect 
is obtained in the sense that a directed current occurs even though there is no net applied 
external biasing force. Unlike many other examples of models of classical ratchets, there is no 
asymmetric potential or asymmetric noise in our model. However asymmetry is incorporated 
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in the modulation of the hopping rates, and this is best seen when we consider the case where 
the modulation is given by vi{t) = sin{ujt—ql). This of course corresponds to a wave traveling 
in a given direction. A non-trivial aspect of the problem studied is the fact that the effect 
goes away as soon as we switch off the hard-core interactions. For non-interacting particles, 
the DC current given by J = (l/LT) Jq dt Yld=i'^iPi ~ '^k+iPi+i, is immediately seen to be 
exactly zero for arbitrary choice of the time-dependent rates. On the other hand, having 
interactions in the system is not a sufficient condition to generate a DC current. For the 
models considered in this paper, the hopping rate is site-wise symmetric. But if the hopping 
rates are symmetric bond-wise, i.e., the hop rate from site / to / + 1 is the same as 

that from / + 1 to /, then the DC current is zero for any choice of phases (pi. To see this, 
consider the density evolution equation obeyed by bond- wise symmetric SEP: 

^ = ui^i,iiPi~i - Pi) + ui,i+iiPi+i - Pi) (44) 

Unlike Eq. ([3]) for site-wise symmetric SEP, p; = p is a solution of the above equation for 
any choice of rates u/. In fact, an inspection of the master equation shows that, even with a 
time-dependent W-matrix, all configurations are equally likely, thus leading to zero current. 
Thus the exclusion process with bond- wise symmetric rates does not give the ratchet effect. 
It is not completely clear as to what are the necessary and sufficient conditions to get a 



directed current 



43|. 



For the model considered here, since the equations for any n-point correlation function 
do not close, it does not seem simple to solve the model exactly. We have therefore studied 
the system analytically using a perturbation theory in the amplitude /i and the inverse 
frequency l/u. In this paper, we have been able to obtain the DC current at order ff by 
solving the evolution equations for density and two point correlation function to order fi. 
This is unlike our earlier solution in 29| where the density was obtained to second order 



in /i. Also, we have been able to obtain results for large driving frequency by solving the 
correlation function alone by such perturbative approaches. Comparing with simulations we 
find that the perturbative results turn out to be quite accurate. 

We now briefly comment on the adiabatic limit, which has been much studied in the 
quantum context. In our case, from our perturbation theory result we see that, over one time 
period of the driving there is a finite particle transport, even in the adiabatic limit. Formally 
we can obtain an exact expression for the net particle transport. For this we start with the 
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master equation dP/dt = W(t)P. Let P^^\t) be the instantaneous equilibrium solution 

satisfying W(t)P(°) = 0. Then, for slow rates to, P{t) will have the form p(°)(t) + uP^^^t) 

where the correction is given by: toP^^^ = dP^'^^/dt . The net particle transported 

across any bond in one time cycle, A/", can then be expressed as: 

"2- ^ ^ dw-L( 



/ dtJ2JiC)P{C,t) = - rfa;^J(C)^^p(°)(C',x) , (45) 
Jo (J •^'^ cc 

where J refers to the current on any given bond. Thus we have a formal expression, for 
the net particle transported, in terms of an integral over an equilibrium average of some 
quantity. However this expression does not appear to have any simple physical interpretation 
and it is not easy to obtain any explicit results, unlike the fast case treated in section fllVp . 
Recently adiabatic pumping phenomena have been studied in the context of geometric phase 



interpretation 



32l |. but the main focus has been on two-state stochastic systems. In this 



case, the current from system to the reservoirs was calculated using full counting statistic 
in the adiabatic or slow driving regime. 

Finally, we point out that an experimental realization of the effect observed in our model 
should be possible in colloidal systems. For instance, consider a colloidal suspension in an 
externally applied laser field. This constitutes a system of diffusive interacting particles in 
an external potential (generated by the laser field) of the form V{x, t) = Vq sin{u!t — qx). This 
system is similar to the model that we have studied. There are some differences, namely, 
in this case because the external field is space dependent, hence the effective hopping rates 
are not symmetric in the forward and backward directions. It would be interesting to study 
this model to see if a current can be generated here, and perhaps one can make detailed 
predictions for experimental observation. 
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